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We have performed the first-principles electronic structure calculation for the novel supercon- 
ductor Ca4Al2 0eFe2As2 which has the smallest a lattice parameter and the largest As height 
from the Fe plane among the Fe-As superconductors. We find that one of the hole-like Fermi 
surfaces is missing around the F point compared to the case of LaFeAsO. Analysis using the 
maximally-localized- Wannier-function technique indicates that the xy orbital becomes more 
localized as the As- Fe-As angle decreases. This induces rearrangement of bands, which results 
in the change of the Fermi-surface topology of Ca4Al2 0eFe2As2 from that of LaFeAsO. The 
strength of electron correlation is also evaluated using the constraint RPA method, and it turns 
out that Ca4Al206Fe2As2 is more correlated than LaFeAsO. 
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The discovery of the new iron-pnictide superconduc- 
tor LaOi-^F^FeAs 1 ) stimulated extensive work to clar- 
ify the basic features of material properties, particularly 
the origin and character of superconductivity and also to 
search for new superconductors. In fact, a variety of re- 
lated superconductors have been synthesized: they have 
the topologically identical iron (Fe)-pnictogen (Pn) layer 
structure combined with various block layers. 2 ~ 4 ' Among 
them, a family of materials which have block layers of the 
perovskite structure are unique in the sense that the dis- 
tance between the neighboring Fe-Pn layers can be sys- 
tematically controlled by changing the thickness of the 
block layer. 4 ) The rich variation of the materials with 
the perovskite structure is also attractive. Quite recently, 
Ogino et al. synthesized (Ca„ + 2(Al,Ti) n 03„_ y )(Fe 2 As2) 
(n=2,3,4) and observed the start of resistivity drop at 
about 39 K for n=4. 5 ' For the synthesis of these materi- 
als at ambient pressure, the presence of Ti turns out to 
be essential not to produce strong internal strain. How- 
ever, subsequent to this work, Shirage et al. have success- 
fully synthesized Ca4Al206-j/Fe2Pn2 with Pn=P and As 
without Ti using a high-pressure synthesis technique. 6 ' 
The superconducting transition temperatures are 17.1 K 
for P and 28.3 K for As. These materials are character- 
ized by the smallest a-latticc parameters and the largest 
pnictogen heights from the Fe plane among the iron- 
pnictide superconductors. It is pointed out that these 
quantities are crucial in determining the transition tem- 
perature. 7 ~ 10 ' 

In the present work, we have performed DFT calcu- 
lations for (Ca4Al 2 06-j/)(Fe2As2) in order to analyze 
the effects of the above unique structural features on 
the electronic structures paying particular attention to 
the energy range near the Fermi level. The absence of 
Ti in this material makes our analysis simpler because 




Fig. 1. Non- magnetic electronic band structure for 
Ca4Al206Fe2As2. The Fermi level corresponds to the en- 
ergy zero. 

d states near the Fermi level comes only from Fe. We 
used a computational code QMAS 11 ' based on the pro- 
jector augmented- wave method, 12 ' which has been ap- 
plied to the study of the ground state properties of 
LaFeAsO 13,14 ' and SrFe 2 As 2 . 15 ' The generalized gradi- 
ent approximation (GGA) 16 ' was adopted for the elec- 
tronic exchange correlation energy. The electronic band 
structure in the vicinity of the Fermi level is analyzed 
in detail using the maximally-localized Wannier func- 
tion technique 17, 18 ' combined with the full-potential lin- 
ear muffin-tin orbital (FP-LMTO) method 19 ' in the local 
density approximation (LDA). 20 ' 

Figures 1 and 2 represent the non-magnetic electronic 
band structure and density of states for the experimen- 
tal structure of Ca4Al206Fc2As2 6 - 1 in GGA obtained by 
QMAS. We found that LDA (using the FP-LMTO code) 
gives almost the same result. As with other iron-pnictide 
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Fig. 2. Electronic density of states (DOS) and partial DOS's for 
Ca4Al20sFe2As2. The Fermi level corresponds to the energy 



superconductors, electronic states near the Fermi level 
derive predominantly from the Fe 3d orbitals though 
there is significant contribution of As 4p states. In Fig. 
3 (a), the Fermi surfaces are shown for the undoped 
case and the doped cases with ± 0.1 electrons /Fe ob- 
tained by the rigid band approximation. There are two 
electron-like surfaces centered at the M point. This is a 
common feature among the iron-pnictide superconduc- 
tors. On the other hand, there are two hole-like surfaces 
around the T point. This is in sharp contrast with other 
iron-pnictide superconductors, where three hole-like sur- 
faces are formed. We will discuss this point later. To 
examine the nesting property, a bare generalized suscep- 
tibility is calculated and plotted in Fig. 3 (b). There is a 
nesting feature from T to M. While the change is small 
with hole doping, it is reduced significantly with elec- 
tron doping, suggesting that the magnetic instability is 
suppressed under electron doping. In the actual samples, 
oxygen content 6 — y~5.6 to 5.8 will produce doping of 
0.2 to 0.4 electrons/Fc. 

In order to analyze the electronic structure in more 
detail, we projected out energy bands near the Fermi en- 
ergy and constructed ten maximally localized Wannicr 
functions. The Wannier functions are characterized by 
the shapes closely related to the Fe 3d orbitals though 
there are significant amplitudes at adjacent As atoms. 
This is consistent with the feature of PDOS's mentioned 
above. The interpolated band structure is then unfolded 
following the procedure of Ref. 21. The result is shown 
and compared with that for LaFcAsO in Figure 4. Here 
the width of the lines represents the weight of each Wan- 
nier function component. In the yz/ zx panel, the weight 
is defined as the average of the weight of d yz and that of 




Fig. 3. (Color) (a) Fermi surfaces for hole-doped (left), undoped 
(middle) and electron-doped (right) Ca4Al20eFc2As2, (b) Bare 
generalized susceptibility for these Fermi surfaces. 



d zx . Our convention for the x and y axes is the same as 
Fig.l in Ref. 21. Namely, we define the x and y for the 
unit cell containing one Fe atom, whereas the X and Y 
refer to the original cell, which is rotated by 45 degrees 
from the x and y axes. 

There are obvious differences in the band structure 
near (0,0) and (tt, it) points between the two com- 
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Fig. 4. (Color online) Unfolded band structure of 
Ca4Al20 6 Fe 2 As2 (left) and LaFeAsO (right). The width 
of the lines represent the weight of each maximally localized 
Wannier function. The Fermi level corresponds to the energy 
zero. The x and y refer to the unit cell containing one iron 
atom in the cell, whereas XY axis is rotated by 45 degrees, 
corresponding to the original cell. 
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pounds in Fig. 4. In order to understand the origins of 
these differences, we calculated the band structure of 
Ca4Ai206Fc2As2 for different values of a while hxing 
both the Fe-As bond length and the distance between 
the As plane and the block layer. The band structure is 
shown in Fig. 5, where the width of the lines represents 
the weight of the xy character. We can see that a does 
affect the band structure significantly. At ci=105 deg., a 
state appears at —0.3 eV at (0,0). This is basically the 
case of Ca4Al 2 06Fc2As2 whose a is 102.13 degrees. As 
a increases, the state moves up and forms the hole-like 
Fermi surface in the case of a=110 deg. With further in- 
crease of a, the band structure becomes similar to that 
of LaFcAsO whose a is 113.55 degrees. The behavior of 
the xy band near (jr, n) with the change of a or the 
pnictogen height from the Fe plane is qualitatively the 
same as what was explained by Kuroki et al. 10 ^ The xy 
band moves up above the Fermi level as a decreases and 
forms a hole pocket near (n, ix), whose presence may be a 
crucial condition for the strong superconductivity of s± 
character. We also found that the distance between FeAs 
layer and block layer is another important parameter. 
With increasing the distance, the band rearrangement 
takes place at a larger value of a. 

As suggested by the above arguments, the xy (X 2 — 
Y 2 ) orbital plays an important role in the band rear- 
rangement. For the analysis of the behavior of this band, 
we list in Table I the onsite energy and the hopping in- 
tegrals between xy (X 2 — Y 2 ) orbitals up to 5th neigh- 
bors, beyond which the hopping integrals are negligibly 
small. We note that the sign of t\ and ti depends on the 
choice of the gauge. We take the opposite sign compared 
to Ref. 22. By cutting off the hopping integrals between 
different kinds of d orbitals, the band dispersion for the 
xy (X 2 — Y 2 ) orbitals is written as 

E xy (k) = e xy + 2ti[cos(fc x ) + cos{k y )} 

+4i2 cos(k x ) cos(fc y ) 

-r-2i 3 [cos(2A: x ) + cos(2A; y )] 

+4<:4[cos(2fc a ;) cos(ky) + cos(k x ) cos(2k y )] 




Fig. 5. (Color online) Unfolded band structure of 
Ca4Al206Fe2As2 for the As-Fe-As angle a =105, 110, 
115 and 120 degrees. The width of the lines represents the 
weight of the xy Wannier function. See the caption of Fig. 4 for 
the convention of the xy axis. 



+4t 5 cos(2k x ) cos(2k y ) , (1) 

where e xy is the onsite energy, and t n is the hopping 
integral for the nth neighbors. Now we can clearly un- 
derstand why the xy band behaves quite differently be- 
tween LaFcAsO and Ca4Al2 0gFc2As2. First of all, the 
onsite energy of Ca4Al2 0eFc2As2 is about 0.4 eV deeper 
than that of LaFcAsO (Table I and Fig. 6) because of 
the change in the crystal field as implied by the dif- 
ference in a. Secondly, the hopping integrals for the 
first three neighbors are quite different between the two 
compounds. Particularly the first neighbor hopping of 
Ca4Al20gFe2As 2 is very small in magnitude and has dif- 
ferent sign from that of LaFcAsO. In Table I, the num- 
bers in the parenthesis for t\ denote the contributions 
from the direct d-d hopping. In the LaFcAsO case, the 
contribution from the indirect Fe-As-Fe hopping (0.394 
eV) overwhelms the direct hopping (—0.243 eV) and 
the net hopping integral is 0.151 eV. However, in the 
Ca4Al20eFc2As2 case, the indirect contribution (0.329 
eV) and the direct contribution (—0.369 eV) nearly can- 
cel each other to produce a small net hopping integral of 
—0.04 eV. The larger direct hopping for Ca4Al2 0eFe2As2 
is due to smaller Fe-Fe distance and the smaller indirect 
hopping is due to smaller a. 

Using Eq. 1 and the parameters in Table I, we found 
that not only the difference in the onsite energy but also 
hopping integrals produce a large difference in the band 
energy at (0,0) between the two compounds. In addition 
to this, the difference in the band energies between (0,0) 
and (7r,7r) is E xy (-n, n) — E xy (0, 0) = —8(^1+2^4). Actual 
values for this quantity are —0.712 eV for LaFeAsO and 
+0.848 eV for Ca4Al20eFe2As2. One may also notice in 
Fig. 4 that the position of the 3z 2 — r 2 branch near (tt, it) 
is quite different between the compounds. For this as- 
pect, we only mention that the trend in the onsite energy 
shown in Fig. 6 is also enhanced by the hopping integrals 
as discussed for the xy (X 2 — Y 2 ) branch. 

Now we discuss briefly other aspects of the band struc- 
ture. As for the x 2 — y 2 (XY) and yz/zx states, the 
hopping integrals are qualitatively the same between 
LaFeAsO and Ca4Al20eFc2As2. However, the variation 
in the onsite energies shown in Fig. 6 may produce qual- 
itatively important effects on the electronic structure 
of these states. In Fig. 6, we first observe the splitting 
between the de(xy, yz, zx) states and the d7(x 2 — y 2 , 
3z 2 — r 2 ) states when the tetragonal distortion is small 
(a~109.4 deg.) in accordance with the ligand field theory 
for the tetrahcdral symmetry. Further splitting among 
states within each of de and d7 as a deviates from 109.4 
deg. is a natural consequence of the distortion in the 
tetrahedral coordination. We also observe a quite inter- 
esting fact that the onsite energy of yz/zx states with 
respect to the Fermi level is nearly invariant for a wide 
range of materials with different as. Other states partic- 
ularly of d7 become deeper and the energy separation 
from yz/zx states becomes larger as a decreases. This 
variation in the onsite energies leads to the increase of 
the weight of de states near the Fermi level as already 
pointed out by our previous papers. 15,22 ) 

Finally we mention the strength of electron corrcla- 



4 J. Phys. Soc. Jpn. 



Letter 



Author Name 



> 
U 

w 



0.4 
0.2 
0.0 
-0.2 
-0.4 
-0.6 



o „ 



— 



- 



"6 u 

3 to 



□ 


x-y (XY) 


• 


yz/zx 


♦ 


3zV 


▲ 





ft i 

is* 



4 
* 



t 

6 



90 100 110 

a (deg.) 



- 

-1 



4 





Ef 



120 



Fig. 6. (Color online) Onsite energy of each Wannicr function for 
eight compounds as a function of the X-Fe-X angle, where X is 
the pnictogen / chalcogen atom. The Fermi level corresponds to 
the energy zero. See the caption of Fig. 4 for the convention of 
the xy axis. 



Table II. Effective Coulomb interaction (U) on the same iron site 
for all the combinations of Wannicr functions in units of eV. The 
definitions of X and Y axes are the same as those in Table I. 
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YZ 


3Z 2 - r 2 


ZX 


X 2 -Y 2 


XY 
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2.21 


2.24 


2.21 


2.53 


YZ 
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2.45 


2.03 


2.07 


3Z 2 - r 2 


2.24 


2.45 


3.46 


2.45 


2.11 


ZX 


2.21 


2.03 


2.45 


2.87 


2.07 


X 2 ~Y 2 


2.53 


2.07 


2.11 


2.07 


2.82 



It would be interesting to see how the pressure changes 
the properties of the compound. Relation with supercon- 
ductivity is also an open question. Magnetic calculations 
including structural optimization are ongoing and will be 
published elsewhere. 
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Table I. Orbital energy and hopping integrals up to fifth neigh- 
bors for xy (X 2 — Y 2 ) state (in eV). The notation with the X 
and Y axes is for easier comparison with the previous works on 
different compounds. 22, 25 ) See the caption of Fig. 4 for the con- 
vention of X and Y axes. For t\, the numbers in parenthesis 
denote the contribution from the direct d-d hopping estimated 
by the dpp model. 
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ti 
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LaFeAsO 


0.189 


0.151 (-0.243) 


0.117 


Ca 4 Al2 06Fe 2 As2 


-0.189 


-0.040 (-0.369) 


0.066 




ts 
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LaFeAsO 


-0.024 


-0.031 


-0.027 


Ca 4 Al20 6 Fe2As2 


-0.012 


-0.033 


-0.029 



tion. The onsite effective Coulomb interaction in the con- 
straint random phase approximation (cRPA) 23,24 ' are 
shown in Table II. The results of other iron-based su- 
perconductors using the same technique can be found 
in Ref. 22,25,26. The average of the diagonal terms is 
3.08 eV, which is larger than that of LaFeAsO (2.53 eV) 
and substantially smaller than that of FeSe (4.24 eV). 22 ) 
The exchange energy (not shown) is smaller than that of 
FeSe, and very close to that of LaFeAsO, with a slightly 
larger value for the xy orbital. The orbital anisotropy is 
weaker than LaFeAsO, reflecting that the xy (X 2 -Y 2 ) 
orbital is less extended. 

In conclusion, we have studied the electronic structure 
of Ca4Al2 0eFe2As2. We found that the topology of the 
Fermi surface is different from that of LaFeAsO. It is re- 
vealed that the As-Fc-As angle has strong impact on the 
band structure in the vicinity of the Fermi level around 
the (0,0) point, specifically through the change in the on- 
site energy level and transfer integral of the xy orbital. 
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